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The works on decoherence due to spin baths usually agree in studying a one-spin system in 
interaction with a large spin bath. In this paper we generalize those models by analyzing a many- 
spin system and by studying decoherence or its suppression in function of the relation between the 
numbers of spins of the system and the bath. This model may help to identify clusters of particles 
unaffected by decoherence, which, as a consequence, can be used to store quantum information. 

I. INTRODUCTION 

o ■ 

(N . 

^ , Decoherence refers to the quantum process that turns a coherent pure state into a decohered mixed state. It is 
^ ■ essential in the account of the emergence of classicality from quantum behavior, since it explains how interference 
' vanishes in an extremely short decoherence time. The orthodox explanation of the phenomenon is given by the 
environment-induced decoherence approach (see Refs. 0, 3) according to which decoherence is a process 
' resulting from the interaction of an open quantum system and its environment. By studying different physical models, 
it is proved that the reduced state pgit) = TrEPsE^J') of the open system rapidly diagonalizes in a well defined pointer 
basis, which identifies the candidates for classical states. 

The environment-induced approach has been extensively applied to many areas of physics —such as atomic physics, 
quantum optics and condensed matter—, and has acquired a great importance in quantum computation, where the 
^ , loss of coherence represents a major difficulty for the implementation of information processing hardware that takes 
' advantage of superpositions. In particular, decoherence resulting from the interaction with nuclear spins is the main 
^ [ obstacle to quantum computations in magnetic systems. This fact has lead to a growing interest in the study of 
O^. decoherence due to spin baths (see Refs. [5] to By beginning from the seminal paper of Zurek ([l[), many works 

have studied the decoherence due to a collection of independent spins. More recently, some papers have directed 
the attention to the interactions between modes within the bath. For instance, by studying a central spin coupled 
I> ■ to a spin-bath, Tessieri and Wilkie ([^) showed that, whereas in the absence of intra-environmental coupling the 
, decoherence of the central spin is fast and irreversible, strong intra-environmental coupling leads to decoherence 
■ suppression. The same model was further analyzed by Dawson et al. (f?!), with the pur pose of relating decoherence 
with the pairwise entanglement between individual bath spins. In turn, Rossini et al. ([iO|) left behind the assumption 
that the central spin is coupled isotropically to all the spins of the bath, and considered the case where the spin system 
interacts with only few spins of the bath. 
. Our analysis can be framed in the context of the above works; it aims at generalizing the paradigmatic spin-bath 
model. In fact, most of the works done so far agree in studying a one-spin system in interaction with a large spin 
bath. The crucial feature of our work is the analysis of a many-spin system, and the study of decoherence or its 
suppression in function of the relation between the numbers of spins of the system and the bath. This generalized 
spin-bath model can also be conceived as a partition of a whole closed system into an open many-spin system and 
its environment. From this perspective, we can study different partitions of the whole system and identify those for 
' which the selected system does not decohere; this might allow to define clusters of particles that can be used to store 
q-bits. 

In order to develop our analysis, we will rely on the general framework for decoherence introduced in Rcf. [l5|, 
where the split of a closed quantum system into an open subsystem and its environment is just conceived as a way of 
selecting a particular space of relevant observables of the whole closed system. Since there are many different spaces 
of relevant observables depending on the observational viewpoint adopted, the same closed system can be decomposed 
in many different ways: each decomposition represents a decision about which degrees of freedom are relevant and 
which can be disregarded in each case. 

On this basis, the paper is organized as follows. In Section II, the standard spin-bath model is presented from the 
general framework perspective: this presentation will allow us to consider two different decompositions, which supply 
the basis for comparing the results obtained for the generalized model in the following sections. In Sections III, IV 
and V, the generalization of the spin-bath model is presented and solved by computer simulations; this task will allow 
us to compare the results obtained for two different ways of splitting the entire closed system into an open system 
and its environment. Finally, in Section VI we introduce our concluding remarks. 
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II. THE SPIN-BATH MODEL 

The spin-bath model is a very simple model that has been exactly solved in previous papers (see Here we 
will recall its main results, obtained from the general framework introduced in Ref. [T5| . in order to compare the 
analogous results to be obtained in the next sections for the generalized model. 



A. Presentation of the model 



Let us consider a closed system U = PUPi U . . . UP/v = -PU (U^^-fj); where (i) P is a spin-1/2 particle represented 
in the Hilbert space tip, and (ii) each Pi is a spin-1/2 particle represented in its Hilbert space Hi- The Hilbert space 
of the composite system U is, then. 



N 



(1) 



In the particle P, the two eigenstates of the spin operator Sp i^ in direction ^ are |-f|-) , such that Sp^ = ^ 
and Sp^^ I-IJ-) = ~5 l-IJ-)- III each particle Pi, the two eigenstates of the corresponding spin operator S^^^ in direction 
it are |ti) , such that S^^^ |ti) = ^ |ti) and \li) = — ^ \li). Therefore, a pure initial state of U reads 



/ N N 

IV'o) = (a 1^) +bm)®{ t.) + /3.I U)) 

\4 = 1 / 



(2) 



where \a\ + \b\ =1 and \ai\ + |/3j| = 1. If the self-Hamiltonians Hp of P and Hi of Pi are taken to be zero, 
and there is no interaction among the Pi, then the total Hamiltonian H of the composite system U is given by the 
interaction between the particle P and each particle Pi (see [l|, [l6|): 



H 



1 ^ 

^(1^) (^i-i^) (^i)®E 



.9. (It.) (td-IW (1.1)® 




(3) 



where = |tj) (tjl + lij) (ijl is the identity operator on the subspace Hj. Under the action of H, the state |?/;( 
evolves into \tp{t)) = a |fr) \£if{t)) + b \£i^{t)) where 



N 

mt)) = m-t)) = (g) (a,e-^'*/2 It.) +/3,e'^-*/2 |;.)) 

i=l 



(4) 



The space O of the observables of the composite system U can be obtained as C = Op (8) {(^fLi^i), where Op is 
the space of the observables of the particle P and Oi is the space of the observables of the particle Pi. Then, an 
observable O € O = % ®'H can be expressed as 



N 



O = Op ® ((g) O.) 



(5) 



where 



Op 



1^) (^1 + %4 1^) (^1 + m (^1 + m 1^) (^1 e 



e^i t.)(t. I + 4ii u){u I + 4ti I + 4? I t.)a. I e 

(0 J«) 



(6) 
(7) 



Since the operators Op and Oi are Hcrmitian, the diagonal components s^^f, sjj,d., e!|.| , e|j^ are real numbers, and the 



(i) 



off-diagonal components are complex numbers satisfying s^^ = s^^, 
observable O in the state \ip{t)) can be computed as 



(i)* 



. Then, the expectation value of the 
(0)v,(,) = (|ap.Si^^ + |6|2s^4)ro(t) + 2Re[a&*s^^ri(i)] (8) 
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where (see [1( 



2=1 

N 

ri(t) = n [l«^l'4^^''* + l/3»^S^"^''* + 2Re(a,/3*6g)] (10) 

2=1 

By contrast to the usual presentations, we will study two different decompositions of the whole closed system U into 
a relevant part and its environment. 

B. The spin-bath model: Decomposition 1 

In the typical presentations of the model, the open system S is the particle P, and the remaining particles Pi play 
the role of the environment E: S = P and E = U^j^P^. Then, the Hilbert space decomposition for this case is 

n = ns®nE^{np)®(^^^'H}j (ii) 



Therefore, the relevant observables Or of the closed system U are those corresponding to the particle P, and they 
are obtained from eqs. (0, ^ and d?]), by making = e|*^^ = 1 and e|j^ = 0: 

Or^Os<»Ie=( J2 I ® ((g)lJ (12) 

The expectation value of these observables in the state |'(/'(i)) is given by 

(Ofl)^(t) = lap + \h\^ HI + 2 Re[a6* s^^ r(t)\ (13) 

where 

r(t) = {£m\£^{t)) = n e-«-* + |/3,p e^^'*) (14) 

i=l 

and, then, 

N 

Vit)? = n(l"*l' + l/^^l' + 2|«.n/32p cos25.t) (15) 

1=1 

This means that, in eq. d!]), ro(i) = 1 and Vx{t) = r{t). 

If we take \oLi\ and as random numbers in the closed interval [0,1], then |r(t)p is an infinite product of 

numbers belonging to the open interval (0,1). As a consequence, limjv^oo r(t) = 0. Therefore, it can be expected 
that, for N finite, r(t) will evolve in time from r(0) = 1 to a very small value (see numerical simulations in [l| and 
0). 

C. The spin-bath model: Decomposition 2 

Although in the usual presentations of the model the open system of interest is P, we can conceive different ways 
of splitting the whole closed system U into an open system S and its environment E. For instance, we can decide to 
observe a particular particle Pj of what was previously considered the environment, and to consider the remaining 
particles as the new environment, in such a way that S ~ Pj and E = P U (U^j^ i-LjPi)- The total Hilbert space of 
the closed composite system U is still given by eq. ([T|), but in this case the corresponding decomposition is 

(16) 
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and the relevant observables Or of the closed system U are those corresponding to the particle Pj : 

Ob ^Os®Ie = Op^®\Ip® (X)^ 1 I (17) 
where (see eq. (O) 

Op, = I t,)(t, \ + ^\ \ + ^^)<t. I + I t,)(i, I (18) 

(?) (?) f?) 

Ip is the identity operator on the subspace Hp, and the coefficients e^^ , , e_|;!|: are now generic. The expectation 
value of the observables Op in the state |'0(i)) is given by 

(Ofl>^(0 = (V'(i)|Ofl,|^(t)) = |a,f e^f + e^^^ + 2 Re (a,/3*e|f e'«*) (19) 

Here there is no need of numerical simulations to see that the third term of eq. (|19p is an oscillating function which, 
as a consequence, has no limit for t — )• oo. This result is not surprising since, in this case, the particle Pj is uncoupled 
to the particles of its environment. 

III. A GENERALIZED SPIN-BATH MODEL: PRESENTATION OF THE MODEL 

Let us consider a closed system U = A\J B where: 

(i) The subsystem A is composed of M spin-1/2 particles Ai, with i = 1, 2, M, each one of them represented in 
its Hilbert space T-La ■ In each A^, the two eigenstates of the spin operator it in direction if are Iffi) and 

Sa„^ \h) = \ \^^) Sa,.^ m = -l m (20) 

AI 

The Hilbert space of A is Ha — ^ ^Ai ■ Then, a pure initial state of A reads 

1=1 

M 
i=l 

(ii) The subsystem B is composed of N spin-1/2 particles Bk, with fc = 1, 2, N, each one of them represented in 
its Hilbert space Hpk- In each B^, the two eigenstates of the spin operator ^ in direction if are |t/c) and 
\ik): 

Sb,,^ \W = \ ItO Sb,^^ \ik) = -\ Mk) (22) 

N 

The Hilbert space oi B y&'Kb — • Then, a pure initial state of B reads 

k=\ 

N 

IV-b) =(8)("fc|tfc)+/3fci;fc)), with \ak\' + \h\' = l (23) 
fc=i 

The Hilbert space of the composite system U = AU B is, then, 

= |^0Ha.^ (24) 



(a) (b) 

FIG. 1: Figure 1: Schema of the interactions among the particles of the open system A (grey circles) and of the open system 
B (white circles): (a) original spin-bath model {M — 1), and (b) generalized spin-bath model {M ^ 1). 



Therefore, from eqs. (|2T|) and (|23l) . a pure initial state of U reads 



iv-o) = iv-a) ® iv-b) = (g) («, 1^,) + h m) U (g) («fe \tk) + Pk life)) 

\4=1 / \fe=l / 



(25) 



As in the original spin-bath model, the self-Hamiltonians and Hsf. are taken to be zero. In turn, there is no 
interaction among the particles Ai nor among the particles B^- As a consequence, the total Hamiltonian H of the 
composite system U is given by 



H = Ha®Hb = 



5fe(ltfe-) (tfel-life) (ifel)® 




(26) 

where — Iftj) (tfjl + l-O-j) (-D-jl is the identity on the subspace Ha^ and Is, = |t/) (t/| + lii) is the identity on 
the subspace "Hs, . Let us notice that the eq. ([3]) of the original model is the particular case of eq. ([26)1 for M — 1. 
This Hamiltonian describes a situation where the particles of A do not interact to each other, the same holds for the 
particles of B, but each particle of A interacts with all the particles of B and vice versa, as shown in Figure 1. 
In eq. (|26p. H is written in its diagonal form; then, the energy eigenvectors are 



i)...|^.) 



h/-i) Y^m) Iti) 
h/-i) \^m) Iti) 



Itfe) ... ItA^-l) M 

ltfe)...|tw-i)M 



h/-i) 14a/) 111)... life) ...i;w-i)M 



(27) 



In turn, the eigenvectors of Ha form a basis of "Ka- In order to simplify the expressions, we will introduce a particular 
arrangement into the set of those vectors, by calling them \Aii: the set {|^i)} is an eigenbasis of Ha with 2^ elements. 
The \Aij will be ordered in terms of the number Z g Nq of particles of A having spin |JJ.). Then, we have that: 



I = corresponds to the unique state with all the particles with spin Ifl"): 

M 

\Ai) = Itt ...,tt) =^ Ha\Ai) = — \Ai) 



(28) 



• 1 = 1 corresponds to the M states with only one particle with spin |-|J-). Since the order of the eigenvectors with 
the same eigenvalue will be irrelevant for the computations, we will order these states in an arbitrary way: 



lA) 

with j 



\i[,i[,...,i[AA:-AA) 

2,3,...,Af + 1 



Ha \A, 



M -2 



lA- 



(29) 



(Af-l)j\/ 



• 1 = 2 corresponds to the ^ 
an arbitrary way: 



states with two particles with spin |JJ-)- Again, we wiU order these states in 



with j = M + 2,M + 3,...,M ■ 



(M - 1) M 
1 + -^^ 



• For the remaining vahies of /, the procedure is analogous. 



Consequently, we have: 



1 eigenvector with eigenvalue 



M 



M eigenvectors with eigenvalue 



M -2 



(30) 



M! . . , . ^ M ~2l 

eigenvectors with eigenvalue 



{M-iy.l 



(31) 



with Z = 0, 1, ...M. Then, it is clear that Ha is degenerate: it has 2*^ eigenvectors but only M different eigenvalues. 
Therefore, a generic state \A) of the system A can be written in the basis as 



|^)-^Q|A)eHA with ^|Q|' = l 

1=1 1=1 

By introducing eq. (15^ into eq. (I^5|) . a pure initial state of the composite system U = AU B reads 
If we group the degrees of freedom of S in a single ket |.B(0)), jV^g) results 



i=l 



The time-evolution of \ip{t)) is ruled by the time-evolution operator U{t) = e = e ^(fiA®HB)t. 



i{HA®HB)t 



mo)) 



If we use Afc to denote the eigenvalue of Ha corresponding to the eigenvector \Ak), then 



where (see eq. (^5]) ) 



|B(t)) = e-'^'^«* |B(0)) 



exp 



AT 



-iAkJ29A\W{tA-\ij){ij\)t 



\m) 



(32) 



(33) 



(34) 



(35) 



(36) 



(37) 
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Since the number of the eigenstates of Ha with the same eigenvalue is given by eqs. (|3T|) . the terms of can be 

arranged as 

m)) = {Ci lA) iSoW» + E 1-^^) i^iW) + E i^a) \B2{t)) 

\\=l J \ X=M+2 

( suaO \ 

+ E Cx\Ax)\Bi{t))\+... + {C2M\A2M)\BMm (38) 

where 

AT 

1^/ W) = (g) (afee^^^'=* Itfe) + /?fce-^^3'=* life)) (39) 
fe=i 

If we compare eq. (p9)) with eq. (|4|), we can see that \£^{t)) and |f||.(i)) are the particular cases of \Bi{t)) for M = 1 
and, then, Z = 0, 1. Let us recall that I is the number of particles of the system A having spin |JJ-)- Then, with M = 1 
and I = 0, \Bi(t)) = \S^[t)), and with Af = 1 and / = 1, \Bi{t)) = \Ei^{t)). 
If we define the function 

rEU(^Oin = o,i,...,M| 

^ ^ 1 otherwise ) ^ ' 



then eq. (|38| can be rewritten as 

Af /(O 

W*))=E E CA|^A)|S;(t)) (41) 

1=0 A=/(i-l) + l 

and the state operator p{t) — \i^{t)) (V'(^)l reads 

/(!) 

A<f /(;') 

-"W^E E ^^AQI^A)|ei(i))(6H0l(^A'| (42) 



i,i'=0 A = /(i-l) + l 
A' = /(i'-l) + l 

An observable OEO—H(E}'Hoi the closed system U — A\J B can be expressed as 

N 



o 



E ^A,A' i^a) (^a'I U (g) (4? It.) (t.i + 4? It.) {u\ + 4? It.) (t.i + 4? it.) a.i) (43) 



, A,A'=0 



Let us notice that eq. ([SJ (a generic observable in the original spin-bath model) is a particular case of this eq. P5|) . 
with only four terms in the first factor. Analogously to that case, the diagonal components sa,a, f-^^l are real 

numbers, and the off-diagonal components are complex numbers satisfying s^.a' — Sy ^tj! — ^it*' Then, the 
expectation value of the observable O in the state p{t) of eq. p2|) can be computed as 

/(!) 

M f(i') 

(0),(t)=Tr(Op(i))= E E ^A,A'rM'(0 (44) 



l,l'=0 A = /(i-l) + l 
A'=/(i'-l) + l 



where 



N 

Ti^i,{t) = Y[ [|a,f e{-f e<*''-»^-.'')5 + e^^'e-^f-'-^,,,')! + 2Re (aj/3*eg'e<9-'+»-'')5)] (45) 



and 



, = (2/ - M) g,, Bx,x' - Ca^^'S^.a (46) 
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Since the exponents in eq. (1451) are of the form gj^i ± gj^v, in some cases they are zero. So, we can write 



M f(M-l) 



M 



/('') 



M f{l) 

(0)pit)=Yl E ^a.a'TmW+E E i?A,A'2Re(T,M-;W)+ E E ^A,A'T,,rW (47) 



1 = A = /(!-l) + l 
A' = /('-l) + l 



1=0 A = /(!-l) + l 

A' = /(M-1-1) + 1 



IJ'=0 A=/(l-l) + l 
Ij^i' A' = /(i-l) + l 
L'^M-L 



where 



M = 



M-2 



if M is even 
if M is odd 



(48) 



T, 



i.i 



N 

w - n [I 



(49) 



TiM-i{t) - n 4t + 1/5/ 4?e-'s„<* + 2Re (a,/3 



(50) 



Let us notice that eqs. ([^^ and ([50]) are analogous to eqs. ^ and pUj) for ro(t) and ri(t), respectively, in the 
original model, with gj i = {21 — M) gj instead of gj. In particular, when M ^ \ and, so, / = 0, 1, then T; i{t) = ro(t) 
andTi,M-i(t) =ri(i).' 

As in the case of the original spin-bath model, here we will consider different meaningful ways of selecting the 
relevant observables. 



IV. GENERALIZED SPIN-BATH MODEL: DECOMPOSITION 1 

A. Selecting the relevant observables 

In this case A is the open system S and B is the environment E. This is a generalization of Decomposition 1 in 
the original spin-bath model. The only difference with respect to that case is that here the system S is composed of 
M >1 particles instead of only one. Then, the decomposition for this case is 



(51) 



Therefore, the relevant observables Or of the closed system U are those corresponding to A, and they are obtained 



from eq. p3)) by making ei*i = ejV = 1, eiV = (compare with eq. ([T2|) in the original spin-bath model): 



N 



Or^Os®Ie=\ E *a.a'|-4a)(-4v| ® (g)I. (52) 

,A,A'=0 / \i=l / 



With this condition, the expectation values of these observables are given by eq. (|47)) . with 

N 

n 

N 



N 



r;j'(t) = n (|a,f e<f^''-f-'')l |/?^.|'e-<9.-.'-9..<')l) (55) 
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4 




/3. 


+ 1/3. 


4 




/3. 



(56) 
(57) 



If we define the functions i?/(t) = \Ti,M-i{t)\'^ and i?;p(t) = |T/,,/(t)|2, they result 

Rlit) = n(i^ 
Rii'it) = n(i' 

i=i 

We can see that |?'(t)p of eq. dTSl) in the original model is the particular case of Ri{t) for M — 1. 

B. Computing the behavior of the relevant expectation values 

The expectation value given by eq. (|47)) has three terms, (0_R)p(t) = S*-^-* + S*^^-* + which can be analyzed 
separately: 

• From eq. ((53)) . the first term reads 

M fil) M f(l) 

^^'^=E E E C^AC^...A^sW(t) (58) 

1=0 A,A'=/(i-l) + l i=0 A,A'=/(i-l) + l 

It is clear that this first term does not evolve with time. 

• The time-dependence of the second term is given by Ti^M^-iit): 

f(i) 

M f{M-l) 

^*'ni)=E E BA,A'2Re(r,,M_Ki)) (59) 

1=0 A = /(i-l) + l 

A' = /(A/-1-1) + 1 



Then, in order to obtain the limit of this term, we have to compute the limit of Ri{t) = |T;_m-;(OP of ^q. 

and \(3j 



2 I 1 2 

As in the case of the original spin-bath model, here we take \aj\ and m J as random numbers in the closed 



interval [0, 1], such that \aj\'' + \l3j\^ = 1. Then 

mayi(\aj\'^ + + 2\aj\'^\Pj\'^cos{2{2l- M)gjt)j (60) 
mm(\aj\* + \pA'^ + 2 \ajf\pA^ cos (2(21- M)gjt)) = f 2 ja^f - (61) 



t 



Therefore, 



I'^il'^ + |/3j I +2|aj|^|/3j| cos {2 {21 — M) gjt) is a random number which, if t 7^ 0, fluctuates 

between 1 and ^2 \aj\^ — ij . Again, when the environment has many particles (that is, when N — > 00), the 

2 I 1 2 2 I 1 2 

statistical value of the cases \aj\ = 1, \f3j\ = 1, \aj\ = and p^ j = tends to zero. In this situation, eq. 



(ISBl) for Ri{t) is an infinite product of numbers belonging to the open interval (0, 1). As a consequence, when 
N ^00, Ri{t) 0. 

• The time-dependence of the third term is given by Ti_ii{t): 

t\i) 

M /(!') 

^^'\t)= J2 E B,,yT,,,it) (62) 



!,!' = A = /(!-l) + l 
i^l' A' = /(i'-l) + l 



with the restrictions on / and I I' and I' ^ M — I. As in the second term, we have to compute the limit 
of Ru'{t) = |T;_i/(t)p of eq. ([57]) and, on the basis of an analogous argument, the result is the same as above: 
when -> 00, Ru'{t) -> 0. 
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If we want now to evaluate the limit of {Ob) p{t) for t — > oo, we have to compute the limits of the second and the 
third terms (since the first term, as we have seen, is time- independent). Here we have to distinguish three cases: 
Af < iV, M > TV and M ~ N. 

Case (a): M < N 

This case is similar to Decomposition 1 in the original spin-bath model, since in both cases M <^ N: the only 
difference is that in the original model M — 1 whereas here M > 1. 

In fact, we have seen that Ti^M-ijt) is analogous to ri(i) in the original model. Moreover, Ti^i'{t) has the same 
functional form as Ti{t). In paper [l6| it is shown that Ti{t) approaches zero for t oo. This means that we can 
infer that Ti^M-iit) and Ti^iit) also approach zero for t — > oo. On the other hand, the terms and TS-^'>{t) 

are sums of less than M terms involving Ti^M-i{t) and Tij'{t). As a consequence, since in this case M is a small 
number, the sum of a small number of terms approaching zero for i — > oo also approaches zero: limf^oo S'^^-'(t) = 
and limt^oo ^'-^\t) = 0. Therefore, 



lim (0;i)p(,) = hm E(i)(i) + + E(3)(t) = j:^^)(t) 



(63) 



In other words, 



M f{l) M f(l) 

1=0 A,A'=/(i-l) + l (=0 A,A'=/(i-l) + l 



(64) 



where is the final diagonal state of U. This result can also be expressed in terms of the reduced density operator 
Pj^ of the system A as: 

hm (0^),(,) = (O^),. = hm (OA)p,(t) = {Oa}p,, (65) 

In the eigenbasis of the Hamiltonian Ha of A, the final reduced density operator is expressed by a 2*^ x 2*^ 
matrix: 



Pa* 














... \ 





Pl=l 








... 








Pl=2 





... 













... 


V o' 











Pi=m) 



(66) 



where P;^q — \Ci\ and each pi is a matrix of dimension ^j^^"^'^,^, x jjjr-f^- This result might seem insufficient 
for decoherence because, since the pi are matrices, seems to be non completely diagonal in the eigenbasis of 
the Hamiltonian Ha- However, we have to recall that all the states \Ai) with same I are degenerate eigenvectors 
corresponding to the same eigenvalue of Ha] then, the basis that diagonalizes (i.e., that diagonalizes all the 
matrices pi) is an eigenbasis of Ha- Summing up, the system S — Aoi M particles in interaction with its environment 
E = B oi N ^ M particles decoheres in the eigenbasis of p^^, which is also an eigenbasis of Ha- 

If we want to compute the time-behavior of {Or} p{t)^ we have to consider that E^-'^^ is a sum of terms of the form 

2 I 1 2 

{B\ \i \aj\ + Bx.y |/3j| ), that is, terms of the expectation value coming from the diagonal part of p{t) in the basis of 
the Hamiltonian H. Therefore, if there is decoherence, the sum Y,"-'^{t) = E^^^ + E^'^\ involving the terms of {Oji) p^t) 
coming from the non-diagonal part of p{t), has to approach zero for t oo. 

In order to show an example of the time-behavior of (0/e)p(f), numerical simulations for E"''(i) have been performed, 
with the following features: 

(i) sa',a = 1 (see eq. (p])). 

(ii) The initial condition for S* = A is selected as (see eq. 



Then, from (i) and (ii), B 



A, A' 




VA, Cx = CI 
(see eq. P5)) ). 



(67) 




0.0002 0.0004 0.0006 0.0008 0.001 

FIG. 2: Figure 2: Evolution of E"'*(t) for N = 10^ and M = 1 (solid line) and M = 10 (dash line), with to = 10~^s. 



(iii) \ai\ is generated by a random-number generator in the interval [0, 1], and |/3j| is obtained as |/3j| = 1 — |ai| . 

(iv) gi — 400 Hz: as explained above, the coupling constant in typical models of spin interaction. 

(v) As in the original model, the time-interval [0,to\ was partitioned into intervals At — to/200, and the function 
S"''(t) was computed at times tk = kAt, with fc = 0, 1, 200. 

(vi) N = 10^, and Af = 1 and M = 10. 
Figure 2 shows the time-evolution of E"'^(t). 

This result shows that, as expected, a small open system S = A oi M particles decoheres in interaction with a large 
environment E = B oi N ^ M particles. 

Case (b): M :$> N 

In this case, where the open system S = A has much more particles than the environment E — B, the argument of 
Case (a) cannot be applied: since now E'-^)(t) and T,^^\t) are no longer sums over a small number of terms, the fact 
that each term approaches zero does not guarantee that the sums also approach zero. In particular, if iV = 1, then 
(see eq. ([55)) ) 

TiM{t) ^ \ai\\<si,i-aij')i + 1/3^2 e"'(si''-9i,'')l (68) 

which clearly has no limit for t — >■ oo. Nevertheless, it might happen that, with N high but M much higher than N, 
each term of the sums approaches zero. So, in order to know the time behavior of (Ofl')p(f), numerical simulations 
for T,"'^{t) have been performed, with the same features as in the previous case, with the exception of condition (vi), 
which was taken as; 



(vi) M = 10^, and TV = 10 and = 100. 




0.002 0.004 0.006 0.008 

FIG. 3: Figure 3: Evolution of E"'^(i) for M = 10^, and N =10 (dash line) and N = 100 (dot line), with to = W'^s. 



Figure 3 shows the time-evolution of in this case. 

This result is also what may be expected: when the open system S = AoiM particles is larger that the environment 
E = B oi N ^ M particles, S does not decohere. 

Case (c): M ~ iV 

In this case, where the numbers of particles of the open system S = A and of the environment = i? do not differ 
in more than one order of magnitude, the time behavior of {Or) p{t) cannot be inferred from the equations. Numerical 
simulations have been performed, with the same features as in Case (b) , with the exception of condition (vi) , which 
was taken as: 

(vi) N =1Q^, and M = 10^ and M = 10^. 

Figure 4 shows the time-evolution of E"'^(t). 

Again, if the environment E = B oi N particles is not large enough when compared with the open system S = A of 
M particles, S does not decohere. Let us notice that, for iV = 10^, the system S = A with M = 10^ does not decohere 
(Figure 4), whereas it does decohere with M = 10 (Figure 2). This shows that, in the case of this decomposition, 
M <^ N means that N is at least two orders of magnitude higher than M. 

Summarizing results 

Up to now, in this Decomposition 1 all the arguments were directed to know whether the system ^ of M particles 
decoheres or not in interaction with the system B of N particles. But, given the symmetry of the whole system, the 
same arguments can be used to decide whether the system B oi N particles decoheres or not in interaction with the 
system ^ of M particles, with analogous results: B decoheres only when M ^ N; ii M <^ N or M N , B does not 
decohere. Therefore, all the results obtained in this section can be summarized as follows: 

(i) li M N, A decoheres and B docs not decohere. 



(ii) E M N, A does not decohere and B decoheres. 




0.001 

FIG. 4: Figure 4: Evolution of Y;"''{t) for iV = 10^ and M = 10^ (dash line) and M = 10^ (solid line), with to = 12.10"^s. 



(iii) If M c:i N, neither A nor B decohere. 



V. GENERALIZED SPIN-BATH MODEL: DECOMPOSITION 2 
A. Selecting the relevant observables 

In this case we decide to observe only one particle of the open system A. This amounts to splitting the closed 
system U into two new subsystems: the open system S is, say, the particle Aj^i with ket Ifl",!!", fl", fl", ff, JJ-), and the 
environment is = {'jfii^A^) \J B = [yjHi^Ai) U (uf^;^Bfc). The decomposition for this case is 



N 



h = Us®He = {nAM)® [y^^ ^-^'j ® V^^""'// ^^^^ 

Therefore, the relevant observables On of the closed system U are those corresponding to the particle Am'- 



Sq,q' |a) (a' 

, Q!,Q'=f|-,JJ. 



N 



(70) 



It is easy to see that the relevant observables selected in this Decomposition 2 form a subspace of the space of the 
relevant observables selected in Decomposition 1: eq. ([70]) can be obtained from eq. (|52p by making s^.a' — 1 for 
A = a' and sa,a' = for A 7^ A' in all the terms of the sum except for the terms corresponding to the particle Am- 
In order to simplify expressions, in this case it is convenient to introduce a new arrangement for the eigenvectors 



of the Hamiltonian Ha, by calling them 



Ai^: the set | -^^)} eigenbasis of Ha with 2*^ elements. The 
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will be ordered by analogy with the binary numbers: 

Ai) = 1^,^,...,^,^,^,^), A2) = \1\A,-AAAA), 



= |^,4,...,4,4,^,4> 



A5) = \i[A,-AAAA), 



A3 
Ae 



A 



According to this arrangement, the 



Ai) with even i have the spin M in the state and the 



the spin M in the state jft)- So, the relevant observables of eq. ([70)) can be rewritten in terms of the 
'^^ = \ Y1 i^n A2\) (A2x + A2X ) (A 



(71) 

Ai ) with odd i have 
Ai) as 



A 



2A-1 



A: 



2\ 



S4JJ, 



A2 



2A-1 / \ A2X-I 



N 



(72) 



B. Computing the behavior of the relevant expectation values 

Here the expectation values of the relevant observables are given by eq. (j47|) . with Ti,ii{t), Tij{t) and Ti^]\j^i{t) 
given by eqs. (|45t . (l49t and (|50)) respectively, but now replacing Bx.y with Bx.x', 



/(O 

A/ f(M-l) 



M 



/(!) 
/('') 



M /(O 

(Ofl)p(*) = E E ^A,A'+E E Kx''2T^^m,M-i{t))+ E E ^a^vTuK*) (73) 



;=o A=/(i-i)+i 

A' = /(i-l) + l 



i=0 A = /(i-l) + l 

A' = /(J\./-i-l) + l 



!,!'=0 A = /(i-l) + l 
Mi' A'=/(i'-l) + l 



where the i?A,A' can be written in the basi; 



is {\Ax)} 



CxCyS^^ if A is an even number and A' = A 

CxCyS-^ij, if A is an even number and A' = A — 1 

^A,A' = { C'aC^ Sjj,^ if A is an odd number and A' = A + 1 

CxCySij.]^ if A is an odd number and A' = A 
otherwise 



According to eq. (|74|) . i?A,A' 7^ only when 



(74) 



A' = A or A' = A ± 1 
Since A = /(/ - 1) + 1 and A' = /(/' - 1) + 1, relations ^ imply that 

I' or I' ^l±l 



(75) 



(76) 



The expectation value given by eq. ([75|) has again three terms, {O) p^t) — S^^^ + S^^^ + which can be analyzed 
separately: 



From eqs. ([7^ and ([75)l . the first term reads 



M f(l) 2""-'' 

^^'^ = E E ^^.^= E {\c2x\'ht + \c2x+i\'hv-)^^^'Ht) 

l=Q A=/(/-l) + l >'=0 



(77) 



Analogously to eq. (|58|) of Decomposition 1, this first term does not evolve with time. 
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The time-dependence of the second term is given by Ti,M-i{t). But with the restrictions of eqs. (|75|) and (|76l) . 
S^^^ has only two terms: 



fC) 

M f{M-l) 

^ E E i3A,A'2Rc(r,^M-iW)= (78) 

i = A = /(i-l) + l 

A' = /(Af-l-l) + l 

= C'^(jw_^_i)+iC'*(Af-i_-^^^2 (^M + 2Re (^Ti^ (79) 
Then, in order to obtain the limit of this term, we have to compute the limit of T m-i m+i (t), which is precisely 

2 ' 2 

the Ti i,{t) of Decomposition 1 in the particular case that I — and V — ^^^^ (see eq. ([55)) ). But, as we have 
seen in Case (a) of Decomposition 1, Ti^i'{t) has the same functional form as ri{t) of the original model (see eq. 
([TU|). which approaches zero for t ^ oo when iV ^ 1. Therefore, for ^ 1, Tm-i m+i (t) also approaches zero 

for t — )• oo, and the same holds for T,^^\t) since it is a sum of two terms containing Tm- i m+i (t). 

• The time-dependence of the third term is given by Ti^ii{t). But with the restrictions of eqs. (|75p and ([76)) . S*^^^ 
results: 

M f(l) 

S<^ni)= E E iBx^\+iTu+i{t) + Bx,x~iTu^iit)) (80) 
Since here I' ^ l±l (see eq. fTG])). in this case Ti^i±i{t) is: 

N 

If we compare this equation with eq. (1141) for r{t) in the original spin-bath model, we can see that 

Tij+i{t)^r{t) and r,,/_i (i) = r*(t) (82) 

Then, 

E(3)(t) = (5+r(t) + 5_r*(t)) (83) 

where 5"+ and S- are constants given by 

M /(O 

^±= E E ^^^^±1 (84) 

, A=/-i A=/(/-i)+i 

2 

On the basis of the simulations of the original model we have seen that, when N ^ 1, r{t) approaches zero for 
t — > oo. Therefore, in this case we can conclude that, when iV ^ 1, E*^'^-'(i) approaches zero for t — > oo. 

Summing up, {Op) p(t) is the sum of three terms: one is time-independent and the other two tend to zero for t — >■ oo. 
In particular, from eq. ([77)) we know that, for iV ^ 1, 

M 1(1) M 2"-^ 

}^{OB)p(t) - E E ^A.A' = E E (|^2a|' ~Sn + l^2A + l|' h^) = {Oe)p, (85) 

/=0 A,A' = /(;-l) + l 1=0 A=0 

where is the final diagonal state of U . Again, this result can also be expressed in terms of the reduced density 
operator Ps = Pam open system S — Am as (see eq. ([65))) 

\im^{OB}p(t) = {Ob}p, = ^hm (OAM)p^,,(t) = {Oam)pa^,, (86) 




0.005 0.01 0.015 0.02 0.025 0.03 

FIG. 5: Figure 5: Evolution of E"'^(t) for M = 10^ and N = 1, with to = S.lQ-^s. 



where the final reduced density operator Pj^^^, in the basis {\-^) , |-IJ.)} reads 

This shows that the open system S = Am, composed of a single particle, decoheres in interaction with its environment 
E of N + M — 1 particles when N ^ 1, independently of the value of M. 

In order to illustrate this conclusion, we have computed Y,"'^{t) = by means of numerical simulations 

with the same features as in Decomposition 1, with the exception of condition (vi), which was taken as: 

Figure 5: (vi) M = 10^ and = 1. 

Figure 6: (vi) M = 10^ and N = 10^ 

Figure 7: (vi) M = 10^ and A'' = 10^. 

Summarizing results 

As we have seen, in this decomposition of the whole closed system, the open system S = Am decoheres when 
N ^ 1, independently of the value of M . But the particle Am was selected as S only for computation simplicity: the 
same argument can be developed for any particle Ai of A. Then, when N ^ 1 and independently of the value of M, 
any particle A^ decoheres in interaction with its environment E of N + M — 1 particles. 

On the other hand, as in Decomposition 1, here the symmetry of the whole system U allows us to draw analogous 
conclusions when the system S is one of the particles of B, say, Sat: S = decoheres when M » 1, independently 
of the value of N. And, on the basis of the same considerations as above, when M ^ 1 and independently of the 
value of A'', any particle Bi decoheres in interaction with its environment E oi N + M — 1 particles. 




VI. CONCLUDING REMARKS 

In this paper we have studied a generahzation of the spin-bath model, where a closed system U is composed by 
two subsystems, U — AU B, with A oi M particles Ai and B oi N particles Bi. We showed how the model behaves 
under different definitions of the system of interest and under different relations between the numbers M and A'^. The 
results so obtained allow us to state the following concluding remarks: 

a) We have seen that, when Af ^ or M ~ N, the subsystem A does not decohere (Decomposition 1 of Section 
IV), but the particles Ai, considered independently, decohere when N ::§> 1 (Decomposition 2 of Section V). This 
means that there are physically meaningful situations, given by M ^ ^ 1 or ~ ^ 1, where all the 
Ai decohere although A does not decohere. In other words, in spite of the fact that certain particles decohere 
and may behave classically, the subsystem composed by all of them retains its quantum nature. We have also 
seen that, by symmetry, all the particles Bi, considered independently, also decohere when M ^ 1. Then, 
when M ^ N ^ 1 oi M c:i N ^ 1, the requirement M ^ 1 holds and we can conclude that not only all the 
Ai, but also all the Bi decohere, although B neither decoheres. So, all the particles of the closed system U — 
{UiAi) U (UjBj) may become classical when considered independently, although the whole system U certainly 
does not decohere and, therefore, retains its quantum character. These results, considered together, are a clear 
manifestation of the fact, already pointed out by Schlosshauer ([i3l), that energy dissipation and decoherence 
are different phenomena: since all the particles of the system U decohere when independently considered, 
decoherence cannot result from the dissipation of energy from the decohered systems to their environments. 

b) The generalized model shows that the split of the entire closed system into an open system and its environment 
amounts to the selection of the observables relevant in each situation. Since there is no privileged or essential 
decomposition, we can select the observables of the subsystem A in the situation in which A does not decohere. 
In this way, it would be possible to use appropriately selected subsystems, unaffected by decoherence, for storing 
quantum information. 

c) The natural further step of generalization will consist in following the ideas of paper , and introducing coupling 
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FIG. 7: Figure 7: Evolution of S"''(t) for M = 10^ and AT = 10^, with to = 4.10-^s. 



internal to the subsystems A or B. For instance, given that the dccoherence of A is increasingly suppressed 
as the number M of its particles increases, it could be expected that such decoherence suppression will also be 
more efBcient as the interactions between the spins of the bath also increase. 
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